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Abstract. We present a novel technique for fitting restframe /-band light curves on a data set of 42 Type la supernovae 
(SNe la). Using the result of the fit, we construct a Hubble diagram with 26 SNe from the subset at 0.01 < z < 0.1. Adding 
two SNe at z ~ 0.5 yields results consistent with a flat A-dominated "concordance universe" (f2 M , Da)=(0.25,0.75). For one of 
these, SN 2000fr, new near infrared data are presented. The high redshift supernova NIR data are also used to test for systematic 
effects in the use of SNe la as distance estimators. A flat, A = 0, universe where the faintness of supernovae at z ~ 0.5 is due to 
grey dust homogeneously distributed in the intergalactic medium is disfavoured based on the high-z Hubble diagram using this 
small data-set. However, the uncertainties are large and no firm conclusion may be drawn. We explore the possibility of setting 
limits on intergalactic dust based on B - I and B - V colour measurements, and conclude that about 20 well measured SNe are 
needed to give statistically significant results. We also show that the high redshift restframe /-band data points are better fit by 
light curve templates that show a prominent second peak, suggesting that they are not intrinsically underluminous. 

Key words. cosmology:observations, supernovae: general 



Send offprint requests to: S. Nobili , serena@lpnhep.in2p3.fr 

* Based in part on observations taken at the European Southern 
Observatory using the ESO Very Large Telescope on Cerro Paranal 
(ESO Program 265.A-5721(B)). Based in part on observations made 
with the NASA/ESA Hubble Space Telescope, obtained at the Space 
Telescope Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under NASA contract 
NAS 5-26555. These observations are associated with program GO- 
8346. Based in part on data collected from the Canada-France-Hawaii 
Telescope Corporation, which is operated by the National Research 
Council of Canada, le Centre National de la Recherche Scientifique 
de France, and the University of Hawaii. Based in part on observa- 
tions taken at the W.M. Keck Observatory, which is operated as a 



1. Introduction 

Observations of Type la supernovae in the restframe fi-band 
at redshifts of z ~ 0.5 and above have shown that they are 
best fit by a cosmological model that includ es a cosmolo gical 
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ported by cross-cutting cosmological results, such as the mea- 
surement of the cosmic microwave backgrou nd anisqtropy, 
which indicates a flat universe (|D e Bernard is et all 12 000: 
IJaffe et allEoOlllSievers et all 120031: ISpergel et alll2003l) : the 
evolution in the number density of X-ray emitting galaxy 
cluste rs JBorgani etail 1200 it iHenrvl EoOll ISchuecker et all 
120031) and galaxy redshift surveys (Efstathio u et all 12002). 
which indicate that Qm ~ 0.3. Taken together, these inde- 
pendent measurements suggest a concordance universe with 
(Q m ,£2a)=(0.25,0.75). However, the SN la Hubble diagram 
remains the most direct approach currently in use for study- 
ing cosmic acceleration, and, thus, possible systematic ef- 
fects affecting the observed brightness of Type la super- 
novae should be carefully consider ed, such as unco rrected 
host galaxy extinction (see e. g. lRowan-Robinsonl (120021) ). dim- 
min g by photon-axion mixing over cosmological distances 
llCsaki et alll2002t|peffaYet et all 1200 lUMortsell et alll2002t 
dstman & Mortsel j,l2004l) and extinction by interga lactic grey 
dust llA guirrel Il999albl: iMortsell & Goobarl l2003h . Some of 
these have alre ady been addres s ed in previous SCP publi- 
cations, see e.g. IPerlmutter et all dl997l fl998l) : LSullivan et alJ 
j2003l) : lKnop et al.N2003l) 

Determining cosmological distances through Type la su- 
pernova fluxes at longer restframe wavelengths offers potential 
advantages, e.g. less extinction by dust along the line of sight, 
either in the host galaxy or in the intergalactic medium. On 
the other hand, the "standard candle" properties at these wave- 
lengths and the possibility of additional systematic effects need 
to be investigated. 

In the restframe /-band, the uncertainties in extinction cor- 
rections are significantly smaller than those in the B-band. For 
example, for Milky-Way type dust (Ry ~ 3) the ratio of ex- 
tinction for the two bands is sizable, Ab/A; ~ 2-3. In gen- 
eral, the extinction corrections become less dependent on our 
knowledge of intrinsic supernova colours and dust properties. 

SNe la /-band light curves typically show a second peak 
15-30 days after the first maximum. It has been suggested that 
the intensity and time-difference between the first and second /- 
band peaks are related to the intrinsic luminosity of the Type la 
SNe, appearing later and more evident for n ormal Type la 
and earlier and fainter for underluminous ones jHamuv et all 
Il996al IWang et all I2005I) . Thus, building /-band light curves 
for Type la supernovae offers the possibility of probing bright- 
ness evolution. 

The scope of this work is to test the feasibility of using 
the restframe /-band for cosmological distance measurements, 
using data available to date, and to assess the importance of ob- 
serving in this wavelength range for future samples of SNe. For 
that purpose, we develop a template fitting technique, which we 
apply to 42 nearby SNe la, to estimate the first (/ max ) and sec- 
ond (/ sec ) /-band light curve peaks. We use the fitted / max of 26 
of these SNe la, which are in the Hubble flow, together with 
two high redshift SNe la: SN 2000fr, at redshift z = 0.543, 
for which new infrared data are presented, a nd SN 1999ff, at 
Z = 0.455, available in literature (Ton rv et all 12003). to build a 
Hubble diagram reaching out to z ~ 0.5. 

The properties of the second peak in the restframe /-band 
light curves are investigated. Furthermore, additional colour in- 



formation is used to test for extinction by non-conventional 
dust fo r the z ~ 0.5 supernovae. In a related work. iRiess et all 
d2000l) used B - I colours of SN 1999Q, in the same red- 
shift range. This SN, however, is not included in our analy- 
sis, since we find inconsistencies with the published data (see 
SectionlO. 



2. l-band light curve fitting 

The second light curve peak seen in /-band for nearby Type la 
SNe varies in strength and position with respect to the primary 
maximum. This complicates the use of a singly parametrised /- 
band template , such as those currently applied in the B- and V- 
band, (see e.g. lPerlmutter et all ill 9971) Jdoldhaber et al.l d200ll) 
for an example of the timescale stretch factor approach), for 
lig ht curve fitting. 

IContardo et ail £2000) proposed a model composed of as 
many as 4 functions for a total of 10 parameters in order to fit 
all UBVRI-bands. Their method used two Gaussian functions 
to fit the two peaks, together with a straight line to fit the late 
time decline and an exponential factor for the pre-max rising 
part of the light curve. In this way, it is possible to describe 
Type la SNe light curves over a wide range of epochs and in 
all optical bands, though, as the authors recognise, it does not 
represent accurately the second peak in the /-band due to the 
influence of the linear decline. However, the main disadvantage 
of their method, for our purpose, is the large number of free 
parameters, which requires very well sampled light curves. 

We have therefore developed a method for fitting /-band 
light curves using five free parameters and one template 1 used 
twice to describe the two peaks. As our goal is only to mea- 
sure the position and amplitude of the two peaks, we limit the 
fit to 40 days after maximum, neglecting the late time decline. 
Our fitting procedure can be summarised as follows: one tem- 
plate is used to fit the time (fi) and the first peak magnitude 
(/1), together with a stretch factor (5/), which is also applied 
to the second template shifted in time to fit the time {to) and 
magnitude of the second peak {It). The underlying function is 

I = hT{s I {t-h)) + hT{s I {t-t 2 )) 

where 7~ is the template. The five parameters fitted are thus: 
If 1 , h, /1 , h, sA, (see Table H). A similar approach is also pro- 
posed by Wan g et alJ (12005) who call it "super-stretch" to em- 
phasise its extension of the stretch approach. 

The use of this fu nction in place of the one described by 
IContardo et alJ d2000t) . reduces the number of free parameters 
by a factor of two. Implicitly, we have thus assumed that the 
rising part of the /-band is the same as in the template used, i.e. 
the B-band. As we will see, this assumption is not always true. 
Note that, unless otherwise specified, the supernova phase al- 
ways refers to the time relative to restframe B-band light curve 
maximum. 



1 The B-band template in Nug ent et alJ 12002) has been used be- 
cause we found that it describes well the data when used with the 
method developed here. We note that there are no physical reasons for 
choosing the B-band over other bands or other kinds of templates. 
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fi time of the peak of the first T template 

Ii peak magnitude of the first T template 

?2 time of the peak of the second T template 

li peak magnitude of the second T template 

si stretch factor of the time axis 

f max time of the first / light curve peak 

^max first / light curve peak magnitude 

f scc time of the second / light curve peak 

/ sct . second / light curve peak magnitude 

Table 1 Summary of the parameters used in this work to de- 
scribe the /-band light curve. The first five parameters are de- 
termined by fitting the data (see text for details). The next four 
parameters are determined from the first set and are the actual 
time and peak values of the light curve. 



2.1. The low-redshift data set 

We applied this method to fit a sample of local SNe la for 
which both B and /-band data are available in the liter ature. 
The SNe prim arily come from th e Calan/Tololo (Hamuv et al., 
Il996ah . CfA jRiess etall Il999h and CfA2 Jjhal l2002h data 
sets. Data from three other well studied individual super - 
novae wer e also included: SN 1 989B dwells et all 1 1994 . 
SN 1994D llRichmond etall 1 19951) and the underluminous su- 
pernova SN 1991bg. We have used two data sets in restframe 
/-band for SN 1991bg, one published bv iFilippenko et ail 
( 1992b) with quite good coverage from about 3 days after B- 
band maximum light to +60 days, and another published by 
iLeibundgut et al.l i ll 9931) with four data points, the first of which 
is at the time of B-band maximum. The agreement between 
the two data sets was assessed by comparing the measurements 
taken at the same date, i.e. JD=2448607, where we found a 
difference of 0.06 mag. We take this as an estimate of the mea- 
surement uncertainty in data of Leibundgut et a i1 lll993l) as no 
uncertainties are reported in that work. 



2.2. Fitting method and results 

Only supernovae with at least 6 /-band data points and time 
coverage constraining both peaks were selected for light curve 
fitting. This resulted in a total of 42 SNe. Table [5]lists the pa- 
rameters resulting from the fitting procedure. Since the domi- 
nant uncertainties are symmetric in units of flux, we performed 
the fit in flux rather than magnitudes. The parameters given in 
Tableware transformed into magnitudes. 

Prior to fitting, all d ata points were ^-corrected to restframe 
/-band as in iKimetalJ l l 19961) and Nugent et al. (2002), assum- 
ing a Bessell /-band filter transmission curve jBessellL Fl990) 
and time information from the available B-band data. A new 
spectroscopic templat e, which is a slightly modified version of 
the template found in Nugent et al. (2002), was built for com- 
puting the ^-corrections. We have preserved the S EP from the 
UV t hrough the Si II 6150A feature, following FNobili et al] 
J2003I) . but red ward of that we have incorporated additional 
spectra from t he Supernova Cosmology Project (SCP) Spring 
1999 search (klderingl 1200(1 iNu gent & Alderingll2000l) to im- 



prove this region as the original template was sparse and re- 
quired a lot of interpolation. A potential source of systematic 
uncertainty in the ^-corrections is due to the wide Ca IR triplet 
absorption feat ure, found to vary c onsiderably a mong Type 
la supernovae dStrolger et all 120021 iNu gent et all |2£)02). We 
have estimated this systematic uncertainty as a function of red- 
shift for 0.01 < z < 0.1 by computing the ^-correction for 
diverse nearby SNe la at different epochs. The dispersion in 
the ^-correction increases with redshift, reaching 0.05 mag. at 
Z = 0.1. We take this as a conservative estimate of the uncer- 
tainty in all ^-corrections in this wavelength region. 

Note that the values of I2 reported in Table [2] are not the 
actual magnitudes of the secondary peak , / sec , but a parameter 
indicating the size of the contribution of the second template to 
the overall I-band light curve. 

In Fig. H a H the fitted light curves are shown. They are 
sorted in chronological order, except for the two very underlu- 
minous supernovae: SN 1991bg and SN 1997cn, displayed at 
the bottom of the figure. As the date of the B-band maximum 
for SN 1997cn is unknown, the origin of the time axis was set 
to the epoch (JD - 2450597.75) when this supernova was first 
observed. Note that the second peak of underluminous super- 
novae is almost completely absent, resulting in a value of I2 ~ 
2.5 to 3 magnitudes fainter than /; . 

Our sample includes SNe that are classified as spectroscop- 
ically peculia r, showing simil a rities w i th the over-luminous 
SN 1 99 IT (IFilippenko et all Il992al iRuiz-Lapuente et all 
ll992HPhillips et allll992fc ~These are SN 1995bd, SN 1997br, 
SN 1998ab, SN 1998es SN 1999aa, SN 1999ac, SN 1999dq 
and S N 1999gn JT.i et alll200ltlHoweiil200lUGaravini et all 
2004). One supernova, SN 1993H, was reported to show 
similarities with t he spec t rum o f the peculiar underluminous 
SN 1986G llHamuv et all Il993l) . However, as we will see in 
this work, we do not find all of these to show peculiarities in 
their /-band light curve shape w hen compared to spectro scop- 
ically normal SNe la. Recently. iRrisciunas et alJ y£)03), built 
the Hubble diagram for SNe la in infrared J, H and /if -bands 
out to z — 0.04, and reported that three spectroscopically pecu- 
liar SNe, SN 1999aa, SN 1999ac and SN 1999aw, do not show 
a behaviour different than that of normal SNe. With the aim 
to assess a greater homogeneity of SNe as standard candles in 
the /-band than in B-band, we choose not to exclude peculiar 
SNe from our sample, and instead monitor possible deviant be- 
haviour of these objects. 

Analysing the results of our fits, we found that Type la SNe 
show a variety of properties for the /-band light curve shape. 
In particular we noticed that the light curve could peak be- 
tween -3 days and +3 days w.r.t. B max , as shown in Fig. |2] 
(left-hand panel). The time of the second peak, f sec (relative 
to B max ), is shown in the right-hand panel. The distribution of 
?max is centred at day -0.3 and has a dispersion of cr = 1.3 
days. f sec is centred at 23.6 with a dispersion of cr — 4.4 days. 
The result show n in Fig. 13 can be compared with the result of 
IContardo et alJ (EoOO) (their Fig. 4). Selecting the subsample 
used by them we obtain a similar distribution, quite flat and 
spread over a broad range, centred around 2 days before B- 
maximum. However, when more SNe are added, we obtain the 
distribution shown in Fig.|2] 
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The fits have reduced^- 2 values (see Tabled) that are gen- 
erally around unity, except for a few cases. 

Approximately half of the fits have reduced x 2 values (see 
Tabled} that are around unity. The other half are either too good 
or too poor, which either suggests that the published uncertain- 
ties are unreliable or that the template is not a good model 
In particular, we note that SN 1994D has a x 2 ld°f ~ 26. 
Although it has been suggested that t he uncertainties fo r this 
supernova may be underestimated (see lKnop et alj j2003)), the 
trend in the residuals shows that this SN is not well described 
by the model, indicating the limitations of the fitting function. 
As in other cases we find a systematic trend, especially in the 
rising part of the light curves for 6 objects, less than half the su- 
pernovae that have pre-maximum data. We note however, that 
a different choice of the template, selected for fitting the pre 
maximum data for these 6 SNe, would fail to fit the rest of the 
sample, which is well fitted by the current template. We in- 
vestigated possible systematic effects in the fitted light curve 
maximum due to this, but found no evidence of a trend in the 
residuals in the Hubble diagram (see SectioiQ for these 6 SNe. 

While the x 2 gives a measurement of the goodness of the 
fit, in the next section we test the robustness and accuracy of the 
parameter estimation in our fitting method, reported in Table|2] 

2.3. Monte-Carlo tests of the fitting method 

Given the heterogeneous origin of the data sample, the quality 
and the sampling of the individual SN light curves vary con- 
siderably. Only a few supernovae have excellent time cover- 
age in the /-band, resulting in a wide range of accuracy in the 
fitted parameters. The robustness of the fitting procedure was 
tested for all circumstances of data quality and time sampling 
in our sample by means of Monte Carlo simulations. We gen- 
erated 1000 sets of simulated light curves for each supernova. 
The synthetic data points had the same time sampling as the 
real light curves and with deviations from the best fit template 
randomly drawn from a Gaussian whose width was set by the 
published uncertainties. The simulated light curves were fitted 
using the same method as the experimental data sets. The dis- 
tribution of the fitted parameters from the simulated data was 
compared with the input data from the fits of the experimen- 
tal data. The mean value in the distribution of each parameter 
generally coincides with that expected, i.e. within one standard 
deviation. There is no evidence for biased fit parameters. This 
lends confidence that the fitting procedure is robust, and given 
the model of the light curve template, will not yield biased es- 
timation of the parameters. 

In two cases, SN 1997br and SN 1998ab, we found that the 
fits to the MC simulations resulted in two solutions, one corre- 
sponding to that found in the fit to the real data and the other 
corresponding to a small fraction (3% and 22% for SN 1997br 
and SN 1998ab respectively) of all simulations. We note, in- 
cidentally, that these SNe are the two with the smallest ratios 
between peak and dip in their light curves. However, the first 
peak is determined by only one and two data points each. A 
close look at the simulated light curves indicates the limited 
number of points constraining the peak is the cause of the rare 
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Fig. 1 /-band light curve fits. On the ordinate is the flux nor- 
malised to the first peak, on the abscissa the restframe time 
since B-band maximum. The dashed line and the dash-dotted 
line represent the two templates used to fit the first and second 
peak respectively. 

failure of the MC simulation. We nevertheless keep these SNe 
in the rest of the analysis, since the parameters and their un- 
certainties estimated from the main distribution agree with the 
results on the real data. 

For the rest of the supernovae, the simulations confirmed 
the expected parameters, giving general confidence in the ro- 
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Table 2 Results of the /-band light curve fit of 42 nearby supernovae: t \ and I\ are the parameters for the time and amplitude fitted 
on the first template, ?2 and h are the parameters for the time and amplitude fitted on the second template, a n d sj is the stretch 
factor. Af is the number of p oi nts used in the fit (do f =N-5). The da ta were taken from: 1 IWells et al ] dl994: 2 iRichmond et ail 
Jl995b : 3 lRiess et alJ Jl999h : 4 lHamuv et ai]dl996al) : 5 ljiialll2002l) : 6 lFilippenko e"ud1ril992bh :lLeibundgut "et alJdl993h . " 



bustness of the procedure and the accuracy of the uncertainties 
on the parameters given in Table |3 

2.4. Intrinsic variations 

We investigated possibl e relations betwe e n /-ba nd and B-band 
parameters. Following iGoldhaber et alJ d200ll) . the time of 
maximum, the stretch factor, sb, and the amplitude of maxi- 
mum, fflg, were determined by fitting a B-band template to the 
published B-band data. A width-luminosity relation was found 
for the first /-band light curve peak. Figure|3]shows the /-band 



absolute magnitude versus the stretch factor in the B-band for 
SNe with zcmb ^ 0.01, where the distance (in Mpc) to each 
SN was calculated from its redshift, assuming a value for the 
Hubble constant, Ho = 72 km s Mpc 1 . The error bars in 
Fig. [3] include an uncertainty of 300 km s _1 on the redshifts to 
account for the peculiar velocities of the host galaxies. The un- 
derluminous supernovae, SN 1997cn and SN 199 lbg, are not 
included in the sample or in any of the analysis presented in this 
section. Corrections for Milky Way and host galaxy extinction 
were also applied, i.e. 

M' mRX - 5 log(// () /72) = / max - Af w - Af st - 25 - 5 \og{d L ) 
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Fig. 1 continued, /-band light curve fits. On the ordinate is the 
the flux normalised to the first peak, on the abscissa the rest- 
frame time since B-band maximum. The dashed line and the 
dash-dotted line represent the two templates used to fit the first 
and second peak respectively. 



Fig. 1 continued, /-band light curve fits. On the ordinate is the 
flux normalised to the first peak, on the abscissa the restframe 
time since B-band maximum. The dashed line and the dash- 
dotted line represent the two templates used to fit the first and 
second peak respectively. 



The host galaxy extinction correction that is applied to most of 
the supernovae is t he weighted average of the three estimates 
given in Table 2 of Phil lips et alJ ( 119991) assuming Rj = 1.82. 
The extinction for the supernovae in the Cf A2 data set was cal- 
culated following the same procedure, using B and V-band pho- 
tometry. At this point, we exclude SN 1995E, SN 1996ai and 



SN 1999cl from the sample as they a re highly reddened (see 
also discussion in Nobil i et al These are not shown 

in any of the plots nor used in any of the analysis that fol- 
low. Two supernovae in the sample, the spectroscopically pe- 
culiar SN 1998es and SN 1999dq, plotted with filled symbols 
in Fig. [5] appear intrinsically redder than average, and become 
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Fig. 1 continued, /-band light curve fits. On the ordinate is the 
flux normalised to the first peak, on the abscissa the restframe 
time since B-band maximum. The dashed line and the dash- 
dotted line represent the two templates used to fit the first and 
second peak respectively. 

~ 2 - 3cr deviant from the average after correction for host 
galaxy extinction. Before introducing light curve shape cor- 
rections, the spread measured in M ! m!lx excluding these two 
SNe, is about 0.24 mag (0.28 mag if they are included). The 
solid line shows the best fit to the data, obtained for a slope 
«/ = 1.18+0.19 and an absolute magnitude for a stretch s B = 1 
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Fig. 1 continued, /-band light curve fits of the underluminous 
supernovae SN 1991bg and SN 1997cn. The dashed line and 
the dash-dotted line represent the two templates used to fit the 
first and second peak respectively. Note that the second peak is 
~ 3 mag fainter than the first peak. 
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Fig. 2 Distribution of the time of /-band maximum referred to 
the time of B-band maximum (left panel) and the distribution 
of the time of second maximum referred to the time of B-band 
maximum (right panel). 



supernova equal to M' max {sB = 1) = -18.89 + 0.03 mag 2 . The 
dispersion, computed as the r.m.s. about the fitted line is 0.17 
± 0.03 mag. A similar correlation was found between the peak 
magnitude and the stretch in the /-band, Si, with a r.m.s. of 
~ 0.19 mag about the best fit line, again excluding SN 1998es 
and SN 1999dq. 

A correlation was found between f sec and the B-band stretch 
factor, as shown in Fig. 0] There are three outliers labelled in 
the figure, SN 1993H , SN 1998es and SN 1999ac, which are 
identified as spectroscopically peculiar supernovae. However, 
other supernovae in our sample that are classified as spectro- 
scopically peculiar behave as "normal" Type la SNe. We note 
that the B-band stretch factor for SN 1999ac is not well defined 
due to an asymmetry of the B-band light curve ( Ph illips et all 
120031) . 

Figure [5] shows a possible correlation between / sec and the 
stretch sb, at least for sb < 0.9, after correcting for the luminos- 
ity distance and for extinction both from host galaxy and Milky 
Way. This correlation, however, disappears for larger values of 
s B - 

All of these correlations, shown in Fig. [3]- [5] were expected 
since it has been suggested that the location and the intensity of 



2 The value fitted for M^ ax depends on the value assumed for the 
Hubble parameter, Ho = 72 km s Mpc -1 . However, its value is not 
used in any of the further analysis presented in this paper. 



Nobili et al.: Restframe /-band Hubble diagram for type la supernovae up to redshift z ~ 0.5 



-20.0 



l 




-18.0 



-17.5 



0.70 



0.80 



0.90 



1 .00 



1.10 



1 .20 



SB 



Fig. 3 /-band absolute magnitude versus stretch in the B-band. 
The best fit gives a 7 = 1.18 + 0.19 and M[ mx (s B = 1) = 
-18.84 + 0.03 mag. The two deviating supernovae, SN 1998es 
(filled diamond) and SN 1999dq (filled circle), were excluded 
from the fit. 



the secondary peak dep ends on the B-b and intrinsic luminosity 
of the supernova jHamuv et all ll996a). 

Figure [6] shows the /-band stretch, sj, plotted versus the 
B-band stretch, sb- We found an interesting linear correlation, 
although some of the supernovae, three of which are spectro- 
scopically peculiar, are more than two standard deviations from 
the fit. The dispersion measured as r.m.s. about the line is 0.08. 

We have investigated the possible existence of further rela- 
tions between the fitted parameters, but find no additional sta- 
tistically significant correlations. 
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Fig. 5 Absolute magnitude of the secondary peak versus the 
stretch in the B-band. The two deviating supernovae in Fig. [3] 
are SN 1998es (filled diamond) and SN 1999dq (filled circle). 
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Fig. 4 Time since B max of the second peak versus the stretch 
in the B-band. The labelled supernovae, classified as spectro- 
scopically peculiar, are excluded from the fit because they are 
outliers. 



Fig. 6 /-band light curve stretch, s/, versus B-band stretch, 
sb- The labelled supernovae are more than two standard devia- 
tion from the correlation shown by the ensemble. SN 1998ab, 
SN 1997br and SN 1998es are classified as spectroscopically 
peculiar. 



3. The /-band Hubble Diagram 

The fitted values of / max were used to build a Hubble diagram 
in the /-band. We select 28 supernovae from the sample con- 
sidered here that have a redshift zcmb ^ 0.01 3 . The maximum 
redshift in this sample is 0.1. 

The width-luminosity relation between I max and the the 
B-band stretch factor was used to correct the peak magni- 
tude, with a a; = 1.18 + 0.19 as measured in the previ- 



3 The lower limit chosen in previous analyses by the SCP is slightly 
higher. However, we include these lower redshift SNe in the sample 
in order to increase the statistical significance. Cutting the Hubble di- 
agram above z = 0.015 would decrease the sample by about 30 %. 
Note, however, that this choice does not significantly affect any of the 
results. 
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Table 3 List of SNe used in the Hubble diagram, mf is the 
peak magnitude corrected for dust extinction and for the width- 
luminosity relation, following Eq. The quoted uncertainties 
do not include the redshift contribution due to peculiar ve- 
locities in the host galaxies, assumed equal to 300 km s . 
Redshifts from Table |2] are here transformed into the CMB 
frame. 



ous section, simi l arly to what is usually done in the fi-band 
JPerlmutter et all [l999). The peak magnitude was also cor- 
rected for Milky Way and host galaxy extinction: 

mf = mi + aj(s B - 1) - A*° " - Af w (1) 

The effective magnitude, mf of the nearby supernovae, listed 
in Table have been used for building the Hubble diagram in 
/-band, shown in Fig. The inner error bars include an un- 
certainty in the redshifts due to peculiar velocities of the host 
galaxies, assumed to be 300 km s _1 . 

The solid line represents the best fit to the data for the con- 
cordance model with fixed Q.m = 0.25 and Qa _ 0.75. The 
singl e fitted parameter, Mi, is defined (as in IPerlmutter et all 
il 19971) 1 to be 

Mi = Mi -5 log H + 25 (2) 

where M/ is the /-band absolute magnitude for a fi-band stretch 
sb = 1 supernova. The value fitted is Mi = -3.19 + 0.03. 
The two redder supernovae, SN 1998es and SN 1999dq, were 



excluded from the fit, and are plotted with different symbols in 
Fig.0 

In order to disentangle the intrinsic dispersion from the sta- 
tistical scatter due to the measurement uncertainties, we sim- 
ulated data sets with a dispersion given by the measurement 
uncertainty only. Since the uncertainty due to peculiar motion 
of the host galaxy is dominant at very low redshift ( ~0.2 mag 
for z =0.01), we limited this calculation to only 15 SNe with 
z > 0.025, which correspond to a peculiar velocity uncertainty 
of the same order as the measurement uncertainties in our sam- 
ple. The average of the r.m.s. measured on each of the sim- 
ulated data sets is geometrically subtracted from the disper- 
sion measured as r.m.s. on the data (0.17 mag), resulting in 
cr — 0.13 mag. We consider this an estimate of the intrinsic dis- 
persion of the stretch corrected /-band l ight curve maximum , 
which agrees with the estimate given bv lHamuv et al.l Jl9 96b) 
using 26 SNe of the Calan/Tololo sample. The estimated in- 
trinsic uncertainty of 0. 13 mag has been added in quadrature to 
the outer error bars of the plotted data. Note that if no correc- 
tion ai(sB - 1) is applied the dispersion in the Hubble diagram 
becomes 0.24 + 0.04 mag, somewhat smaller than the cor- 
responding dispersion measured in the "uncorrected" fi-band 
Hubble diagram. Moreover, we computed the dispersion in the 
Hubble diagram for the three data sets separately, and no sta- 
tistically significant differences were found. 

4. High redshift supernovae 

Next, we explore the possibility of extending the Hubble dia- 
gram to higher redshifts, where the effects of the energy density 
components of the universe are, in principle, measurable. The 
restframe /-band data available to date for this purpose are un- 
fortunately very limited. They consist of only three supernovae 
(SN 1999Q, SN 1999ff and SN 2000fr) at redshift z ~ 0.5 ob- 
served in the near infrared (NIR) J-band collected during three 
different campaigns conducted using different facilities and by 
two different teams. Keeping all of these possible sources of 
systematic errors in mind, we include two of these supernovae 
in the /-band Hubble diagram. 

4.1. SN2000fr 

SN 2000fr was discovered by the Supernova Cosmology 
Project (SCP) during a search for Type la supernovae at 
redshift z ~ 1 conducted in the /-band with the CFH12k 
camera on the Canada-Fr ance-Hawaii Telescope (CFHT) 
dSchahmaneche et"all Eo01). The depth of the search allowed 
us to discover this z ~ 0.5 supernova during its rise, about 1 1 
rest-frame days before maximum B-band light. 

The supernova type was confirmed with spectra taken at the 
Keck II telescope and th e VLT, showing that it was a no rmal 
Type l a at z = 0.543 ( see Lidman et al.1 ( EooH) : iGaravini et alJ 
(2005J) for an extensive analysis of the spectrum). This super- 
nova was followed in the restframe B, V and / bands involving 
both ground and space based facilities. Approximately one year 
later, when SN 2000fr had faded sufficiently, infrared and opti- 
cal images of the hos t galax y were obtained. The optical light 
curves in Kno o et alJ (|2Q03) were re-fitted using the improved 
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Fig. 7 Effective /-band maximum vs redshift for the nearby supernovae of the Calan/Tololo, CfA and Cf A2 sample. The data have 
been corrected for the stretch-luminosity relation and for Milky Way and host galaxy extinction. The r.m.s. along the concordance 
model line is cr - 0.17 + 0.04 mag. Subtracting the contribution of the average uncertainty, results in 0.13 mag estimated intrinsic 
dispersion (see text for details). SN 1998es and SN 1999dq were excluded from the fit (see text) and are plotted with open 
diamonds. The inset plots the histogram of the residuals. 



spectral templates for computing ^-corrections. We found a B- 
band stretch factor of s# = 1.034 + 0.013 and a time of Z? max , 
Wx = MJD 51685.6. Restframe B - V measurements at the 
time of fi m ax indicate that SN 2000fr did not suffer from red- 
dening due to dust in the host galaxy (see Section[6]for a more 
extensive discussion). The adopted M ilky Way reddening is 
E(B - V) = 0.030 mag ISchlegel et allfl998i). 

The near-infrared data were collected with ISAAC at the 
VLT. They consist of 7 s -band observations during three epochs 
and a final image of the host galaxy without the SN (see Table 
0}. Each data point is composed of a series of 20 to 60 images 
with random offsets between exposures. Figure[8]shows a com- 
parison between the Persson J filter and the narrower ISAAC J s 
filter used for the observations, together with the atmospheric 
transmission, and the spectral template at maximum. 



The advantage in using the narrower J s filter is that the 
transmission of the filter is not determined by the region of 
strong atmospheric absorption between 13500 and 15000 A. 
Consequently, the zero-point is significantly more stable than 
that of standard J. This was very useful, because all the ISAAC 
data were taken in queue mode, where typically only one or two 
standard stars, chosen from the list of Persson e tai]dl998l) . are 
observed during a night. All data, except the reference images, 
were taken during photometric nights and the difference in the 
zero-points from one night to the next was less than 0.01 mag- 
nitudes. 

The data were reduced using both internally developed rou- 
tines and the XDIMSUM package in IRAF 4 . The differences 
between the two analyses are within the quoted uncertainties. 

4 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
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Fig. 8 Comparison between the Persson / filter (dashed-dotted 
line), the J s filter at ISAAC used for the observations of 
SN 2000fr (dashed line) and the /-band red-shifted to z=0.543 
(heavy solid line). The atmospheric transmission is also plotted 
(dotted line). The spectral template at day is on an arbitrary 
flux scale for readability purpose (solid line). 



MJD 


Epoch 


/,(mag) 


/(mag) 


51685.06 


-0.33 


22.50 ± 0.09 


23.52 ±0.10 


51709.02 


15.20 


23.57 ± 0.22 


24.52 ± 0.23 


51731.96 


30.07 


23.14 ±0.15 


23.99 ±0.16 



Table 4 Summary of / s -band data for SN 2000fr. The quoted 
errors are due to statistical Poisson noise and the uncertainty in 
the ZP (contributing 0.01 mag). Epochs are in restframe days 
relative to the date of B-band maximum. The restframe /-band 
is obtained through cross-filter ^-correction from the observed 
/ s -band to Bessel /-band. The uncertainties also include the 
contribution from ^-corrections, estimated to be 0.05 mag at all 
epochs considered. 



The supernova images were aligned with the host galaxy im- 
ages and the flux scaled to the one with best se eing, using th e 
field stars before performing PSF photometry (Fabbro, 2001). 
The results are presented in Table |4] The stated uncertainties 
include the statistical Poisson noise and the uncertainty on the 
estimate of the zero point, added in quadrature. 

The / s -band magnitude takes into account a colour term 
which arises from the difference between the / filter of the stan- 
dard star system and the J s filter used in ISAAC. This correc- 
tion was small, ~ 0.012 mag. 

The cross-filter ^-correction, Kjj s , to convert fro m 7, -band 
to rest -frame /-band, has been calculated following Kim et all 
(1996) using the spectral templates improved for this work. 
The ^-correction includes a term to account for the appropriate 
transformation between IR and optical photometric systems, 
equal to (I — J) = 0.03(+0.02), determined by using the Vega 



magnitudes in / and / jBessell et allll998HCohen et alll2003h . 
We conservatively assume 0.05 mag total uncertainty in the k- 
corrections (see Section l2~2l . 

4.2. SN 1999ff 

SN 1999ff was discovered by the High-Z Supernova Search 
Team (HZSST) during a search cond ucted at C FHT using the 
CFH12k camera in the /-band (T onrv et"aill2003l) . 5 The super- 
nova was confirmed spectroscopically as a Type la at redshift 
z = 0.455. The adopted Milky W ay reddening is E(B - V) = 
0.025 mag iSchlerel et allfl998i). 

/-band observations, cor responding to restframe /-band, 
reported in To nrv et al J ( 120031) . were taken at Keck using NIRC 
at two epochs only. The /-band filter that was used for these 
observations is very similar to the ISAAC J s , shown in Fig. (HI 
We have used the published photometry, and, for consistency 
with the treatment of both the low redshift supernovae and 
SN 2000fr, we computed the ^-corrections using the improved 
spectral templ ates. We foun d differences with the results pub- 
lished in Tonr v et ail d2003h . due to the use of an incorrect filter 
in the originally published results. (However, the ^-corrections 
calculated as part of the MLCS distance fits to this object were 
done with the correct filter (Brian Schmidt, private communi- 
cation)). The /-band magnitudes were also corrected for the 
offset found between the optical and IR systems, as explained 
in the previous section. The restframe /-band magnitudes ob- 
tained this way are reported in Table [5] The published optical 
B-band data were used to fit the restframe B-band light curve 
using the stretch method. Our t ime of maximum was within 
1 day of the iTonrv eta fl d2003l) value, with a best fit for the 
stretch s B = 0.80 + 0.05. 



MJD 


Epoch 


I(mag) 


51501.29 


5.01 


23.57 ±0.11 


51526.31 


22.21 


24.06 ± 0.24 



Table 5 Summary of IR data for SN 1999ff. Epochs are in 
restframe days relative to the date of the B-band maximum 
(fmax = MJD5 1494.8); restframe /-band magnitudes are com- 
puted applyin g ^-correction s to the observed /-band data pub- 
lished in Tonr v et alJ (|2003). The uncertainties also include the 
contribution from ^-corrections, estimated to be 0.05 mag at all 
epochs considered. 



4.3. SN 1999Q 

SN 1999Q was discovered by the HZSST using the CTIO 4m 
Blanco Telescope and was spectrally con firmed to be a Type la 
SN at z = 0.46 (Gamavidi£t^l|l999). The adopted Milky 



for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



5 Another supernova, SN 1999fn, was followed in /-band by the 
HZSST during the same search. However since it was found in a 
highly extincted Galactic field, E(B-V)=0.32 mag, and since it was 
strongly contaminated by the host galaxy, we did not include it in our 
analysis. 
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Way reddening is E(B-V) = 0.021 magnitudes dSchlegel et all 

SN 1999Q was observed in the /-band over five epochs, 
the first with Sofl on the ESO NT T and the following four with 
NIRC at the Keck Telescope jRiess et all 120001) . We recom- 
puted the ^-corrections using our new spectral template (as we 
did for SNe 2000fr and 1999ff) and we find a difference of up 
to 0.15 magnitudes betwee n our ^-corrections and those pub- 
lished in lRiess et allfcOOOl) . 

A fit to the published restframe /-band data of SN 1999Q 
shows that it is a 4 standard deviation outlier in the /-band 
Hubble diagram. In order to investigate its faintness, we re- 
analysed the publicly available Sofl data and found 7=22.63 
+ 0.15 mag, which is sig nificantly br ighter than the published 
value, 23.00 + 0.14 mag dRiess et allboOOh . Due to this large 
discrepancy, we decided to not include this SN in the rest of the 
analysis. 

4.4. Light Curve fits for the high redshift supernovae. 

The /-band light curves of the high redshift supernovae are not 
as well sampled in time as the low redshift sample analysed in 
Section|2 There are only few data points for each SN, making 
it impossible to perform the 5 parameter fit. Thus, we used the 
results of the fit of the local sample of supernovae to build a set 
of 42 /-band templates, which in turn have been used to fit the 
high redshift SN light curves. 

The best fit light curve for each of the 42 supernovae in 
our low-redshift sample can be viewed as defining an /-band 
template. The high redshift supernovae are fit to each template 
with a single free parameter, / max , the absolute normalisation 
of the template. The time of B max is obtained from the lit- 
erature (SN 1999ff) or from our own B-band light curve fits 
(SN 2000fr). The best-fitting low-redshift /-band template fixes 
the date of the I-band maximum relative to the date of the B- 
band maximum. A x 1 comparison was used to choose the best 
low redshift template. Fi g s . l9l and 1 1 01 sho w the comparison of 
the data with the best fit template for each of the supernovae. 
Table [6] gives the results of the fit together with redshift, the 
number of data points, the template giving the best fit and the 
X 1 - As there are only a few data points for each SN, the^ 2 pa- 
rameter has little significance for estimating the goodness of 
the fits. Thus, to estimate the possible systematic error in the 
measured peak magnitude from the selection of the light curve 
template, we computed the r.m.s. of the fitted /„ MA of all the 
light curve templates satisfying \ 2 — xlm, + 3- This possible 
systematic uncertainty is reported also in Table[6] For both SNe 
this is quite small, and compatible with the scatter due to the 
statistical uncertainties, thus, it is a conservative estimate. 



4.5. Monte-Carlo test of the fitting method 

A Monte-Carlo simulation was run in order to test the robust- 
ness of the fitting method applied to the high redshift SNe. The 
measurement uncertainties were used to generate a set of 1000 
SNe, with data points randomly distributed around the real data 
and at the same epochs as the data. All the simulated data sets 



;| 0.6 

a 




10 20 30 

Days since Bmax 

Fig. 9 /-band fit for SN 2000fr. Out of 42 /-band templates, the 
best fit was obtained with the template of SN 1992bc. The fit 
was performed with only one free parameter, the peak magni- 
tude, / max = 23.48 + 0.08 mag. Supplemental data from the 
B-band (not shown) is used to fix the date of B-maximum. 




10 20 30 

Days since Bmax 

Fig. 10 /-band fit for SN 1999ff. Out of 42 /-band templates, 
the best fit was obtained with the template of SN 1996bl. The 
fit was performed with only one free parameter, the peak mag- 
nitude, / max = 23.55 + 0.10 mag. Supplemental data from the 
B-band (not shown) is used to fix the date of B-maximum. 



were in turn fitted with the 42 templates and the one giving 
the minimum x 1 was selected for each simulation. The distri- 
bution of the fitted parameters in each of the simulated data 
sets around the true values, fitted on the experimental data, was 
studied to check for systematic uncertainty in the fitting pro- 
cedure. This was found to be robust, always selecting the same 
template as the one giving the best fit for both SNe. No bias was 
found, therefore confirming the peak magnitude fitted with this 
method. The uncertainty in / max reported in Table [6] was con- 
sistent with the dispersion in the distribution of / max measured 
from the simulations. 



5. The /-band Hubble diagram up to z ~ 0.5 

The /-band peak magnitudes of the high redshift supernovae 
reported in Table [6] were corrected for Milky Way extinction. 
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SN 


■7 


SB 


n 








template 


1 

X~ 




SN 2000fr 
SN 1999ff 


0.543 
0.455 


1.034 ± 0.011 
0.80 ± 0.05 


3 
2 


23.48 
23.55 


+ 0.08 
+ 0.10 


±0.04 
±0.08 


SN 1992bc 
SN 1996b! 


1.04 
0.05 


0.027 
0.022 



Table 6 List of the high redshift Type la SNe used in this work. Columns are: IAU name, redshift, number of data points used 
in the fit, magnitude of the peak resulted from the fit (both statistical and systematic uncertainties are given) before Milky Way 
extinction correction, best fit template, x 1 of the fit, Milky Way extinction in the /-band. 
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Fig. 1 1 Effective /-band maximum versus redshift for the nearby supernovae of the Calan/Tololo, CfA and CfA2 sample, together 
with two supernovae at redshift ~ 0.5 for case c (top panel), residuals to the (Qm, Oa)=(0.25,0.75) model for case c (middle 
panel) and case a (bottom panel). SN 1998es and SN 1999dq were excluded from the fit (see text) and are plotted with open 
diamonds. Only the low-redshift SNe have been extinction corrected. The host galaxy extinction correction would make the two 
high-redshift data points closer to each other but with larger error bars. 



Note that both SN 1999ff and SN 2000fr have been reported 
not to suffer from extinction from their host galaxies. 

The Hubble diagram has been built both with and without 
width-luminosity correction (case c and case a respectively), 
where the systematic uncertainties on the peak magnitudes of 
the distant supernovae, listed in Table |6j are added in quadra- 
ture to the statistical uncertainties. Cases b and d are like a 
and c but neglect the systematic uncertainties from Table [6] 



Figure^2 snows the extended Hubble diagram (case c), where 
an intrinsic uncertainty of 0. 13 mag has been added in quadra- 
ture to the measurement errors of the plotted data. The solid 
line represent the best fit to the nearby data for the concor- 
dance model Q.m = 0.25 and Qa = 0.75. Also plotted is the 
model for Q.m = 1 and Qa = (dashed line), and a flat, A = 
universe in the presence of a homogeneous population of large 
dust grains in the intergalactic (IG) medium able to explain the 
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observed dimming of Type la SNe at z ~ 0.5 in the B-band 
(dotted line) ( Aguirre, 1999a b). The bottom panel shows the 
residuals obtained for case a. Table[7]hsts the^ 2 values for the 
high redshift SNe for each of the models. The £2 A -dominated 
cosmology is formally favored over the other two models at 
the >2cr level. However, two high redshift supernovae obvi- 
ously do not provide the full gaussian distribution that would 
confirm this result. 



together (see Table II Oi l. The correlati ons between SN colours 
at different epochs found in (Nobili et al., 2003) were taken 
into account. However, we note that, although this correlation 
should be taken into account in the calculations, neglecting it 
would not change significantly the conclusions of the analy- 
sis. Although individual supernovae give x 2 values that would 
seem to distinguish between the models, the combined results 
disfavour such conclusions. 



(fl M , Ha) 


Xa 


xl 




Xd 


(0.25,0.75) 


2.18 


2.36 


1.16 


1.44 


(1,0) 


7.56 


8.44 


15.56 


18.30 


(l,0) fa 


3.52 


3.96 


5.36 


6.48 



Table 7 x 2 (f° r 2 dof) of each model to the high redshift data, 
without stretch correction and with systematic uncertainties 
added in quadrature (x 2 ), neglecting the systematic uncertain- 
ties (xl), with stretch correction and adding the systematic un- 
certainties in quadrature (x 2 ) or neglecting them (x^)- 



Epoch 


B-I 


SN 2000fr 


-0.32 


-0.51 ±0.12 


14.70 


-0.50 ± 0.24 


29.08 


1.42 ±0.17 


SN 1999ff 


5.59 


-0.18 ±0.13 


27.20 


1.38 ± 0.24 



Table 8 Restframe B-I colours in magnitudes for the two high 
redshift SNe. The Epoch is in restframe days relative to the B- 
band maximum, divided by the B-band stretch. 



Systematic uncertainties in the method used here also can- 
not be extensively explored with only two supernovae. Some 
uncertainties are specific to the sample considered here. The 
different fitting methods applied to the restframe /-band light 
curve for the low and high redshift samples can be easily over- 
come if distant supernovae are followed at NIR wavelengths 
with better time coverage. Both the low and high redshift sam- 
ples used in this analysis are rather heterogeneous, as they were 
collected from different data sets. Future data sets collected 
with a single instrument would naturally solve this problem. 

6. SN la colours and intergalactic dust 

Multi-colour photometry allows one to search for non-standard 
dust having only a weak wavelength dependence, such as a 
homogeneous population of large grain dust, as proposed by 
lAguirrel(ll999alhl) . 

If we assume that grey dust is responsible for the dim- 
ming of SNe la in the B-band at z ~ 0.5, we can calculate the 
expected extinction in other filter s and c ompute the resulting 
colours. Following iGoobar et alJ d2 002al). w e use the SNOC 
Monte-Carlo package dGoobar et alT B002b) for two cases of 
the total to selective extinction ratio - Ry — 4.5 and 9.5. We 
assume that the dust is evenly distributed between us and the 
SNe in question and we assume a flat cosmological model with 
a zero cosmological constant. 

The measured B-I and B —V colours of SN 1999ff and 
SN 2000fr, corrected only for Milky- Way extinction, are pre- 
sented in Tables[H]and[9]and plotted in figure^l The error bars 
include the contribution of the intrinsic colour dispersion. The 
expected evolution in the B-I and B — V colours of an average 
SNe la in the concordance model and in the two models with 
grey dust and without a cosmological constant at z — 0.5 are 
also shown. 

The x 2 has been computed for both B — V and B-I evo- 
lution for SN 1999ff and SN 2000fr, and for both supernovae 



Epoch 


B - V 


SN 2000fr 


-7.97 


-0.06 ± 0.05 a 


-3.51 


-0.14±0.05 a 


4.60 


-0.12 ±0.05 


12.93 


0.24 ± 0.08 


20.31 


0.61 ± 0.07 


30.22 


0.99 ± 0.09 


SN 1999ff 


-7.99 


0.03 ± 0.08 a 


1.91 


-0.02 ± 0.09 


1.98 


0.10 ±0.12 


2.91 


0.23 ± 0.12 


19.55 


0.71 ±0.09 


28.75 


1.22 ±0.20 



Table 9 Restframe B - V colours in magnitudes for SN 1999ff 
and SN 2000fr. The Epoch is in rest frame days relative to the 
B-band maximum, divided by the B-band stretch. 



"The data are not inc luded in the a nalysis because they are out of the 
range in which lNobili et alJ 120031) studied colour correlations. 



To make our test for grey dust more effective, a different 
approach was followed. The method of least squares has been 
used to combine colour measurements along time for each su- 
pernova (see Cowan, 1998, p. 106 for details). The residuals be- 
tween the data and the models are averaged with a weight that 
is determined from the covariance matrix. In the following, we 
refer to E(X - Y) to describe the colour excess of any super- 
nova with resp ect to the average X - Y colour of nearby SNe la, 
as derived in (iNobili et all E003). First we applied this method 
to all local supernovae and used the results to establish the ex- 
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Fig. 12 The evolution in the colour of SN 1999ff (squares) and 
SN 2000fr (triangles), B — V (top panel) and B - I (bottom 
panel), compared to the colour evolution of the average SNe la 
in a A dominated universe (solid line) and a Q.m = 1, = 
universe with presence of grey dust with Ry=A.5 (dashed line) 
and Ry-9.5 (dotted line) forz = 0.5. 



pected distribution in the E{B-I)\s E(B- V) plane, as showed 
in Fig. El 

As the high redshift SNe were not corrected for host galaxy 
extinction, we computed the colour distribution of nearby 
SNe la for two cases: the left-hand panels represent the dis- 
tribution of colour excess of 27 nearby SNe not corrected (top 
panel) and corrected (bottom panel) for host galaxy extinction. 
Spectroscopically peculiar SNe have been excluded from the 
analysis. The projection of the ellipses on each colour axis is 
the estimated standard deviation in that colour and the incli- 
nation is defined by the linear Pearson correlation coefficient 
computed on the same data sample. The solid contours repre- 
sent 68.3%, 95.5% and 99.7% probability. 

The right-hand panels in Fig. II 31 show the combined val- 
ues of colour excess, [E(B - V),E{B - I)], for the high red- 
shift supernovae: [0.12 + 0.09,0.25 + 0.17] for SN 1999ff and 
[-0.1 1 +0.08, -0.10+0.18] for SN 2000fr. These are compared 



to the local supernova distribution (solid lines), that represent 
the distribution expected in the absence of IG dust. Also plot- 
ted is the 68.3% level of the expected distribution in presence 
of grey dust with Ry = 9.5, represented by the ellipse (dashed 
line) that is displaced by (0.06,0.19) from the no-dust model. 
Only the case of Ry = 9.5 has been plotted for readability rea- 
sons, given the small difference between the two dust models. 
Note that this is the closer to the no-dust model. The ellipse 
corresponding to Ry = 4.5 would be displaced by (0.03,0.04), 
respectively in E(B-V) and E(B-I), from the Ry = 9.5 model. 

We computed the y 2 of the high redshift data for all three 
models, for two cases: in the first case, the nearby SNe la are 
corrected for extinction by dust in the host galaxy, and, in the 
second case, they are not (bottom and top panels of Fig. II 31 . 
For each model, we sum the y 1 contribution from all SNe, tak- 
ing into account the correlation found between E(B - V) and 
E(B - I) in the nearby sample. In the first case, the reduced^ 2 
(for 4 degrees of freedom) are 0.63, 1.20, and 0.97 for the no- 
dust, IG dust with R v - 4.5 and IG dust with R v = 9.5 models 
respectively. In the second case, the reduced y 1 are 0.62, 1 .55 
and 1.32 respectively. We note that both the intrinsic disper- 
sion in the colours of the nearby data, and the uncertainties in 
the colours of the high redshift SNe, have been taken into ac- 
count in computing the^ 2 . The statistical significance of these 
results is very limited, and should only be taken as an exam- 
ple of the method developed here. Moreover, the possibility for 
this analysis to be affected by systematic effects is not negli- 
gible. Increasing the sample and the time sampling for each 
object would allow us not only to improve the significance of 
our statistic, but it will also be a means to identify and quantify 
systematic effects involved. 

A Monte Carlo simulation was used to estimate the mini- 
mum sample size needed to test for the presence of homoge- 
neously distributed grey dust in the IGM. SNe colours were 
generated following the binormal distribution defined by the 
nearby SNe la sample. Under the assumption that the system- 
atic effects are negligible and an average measurement uncer- 
tainty of 0.05 mag in both E(B - V) and E(B-I), we found that 
a sample of at least 20 SNe would be needed to exclude the IG 
dust model with R v = 9.5 at the 95% C.L. Note that the average 
measurement uncertainty of 0.05 mag can be achieved with dif- 
ferent strategies. Currently, the uncertainties on the individual 
measurements give the main contribution to the colour uncer- 
tainties. A good sampling would allow us to better identify and 
quantify currently unidentified systematic effects which may 
possibly be affecting the current analysis. 

7. Test for SN brightness evolution 

Evolution of the properties of the progenitors of SNe la with 
redshift has often been proposed as an alternative explanation 
for the observed dimming of distant SNe. This is based on the 
assumption that older galaxies show different composition dis- 
tribution than younger ones, e.g. an increased average metal- 
licity, resulting in different environmental conditions for the 
exploding star. A simple way to test for evolution is to com- 
pare properties of nearby SNe with distant ones. This will not 
prove that there is no evolution, but it will exclude it on a 
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Fig. 13 Left-hand panels: distribution of combined colour measurements of nearby SNe la in the E(B - 1) vs E(B - V) plane, not 
corrected (top panel) and corrected (bottom panel) for the host galaxy extinction. The solid contours incorporate 68.3%, 95.5% 
and 99.7% of the sample. Right-hand panels: SN 1999ff and SN 2000fr compared to the distribution defined by nearby SNe in 
the case of no IG dust (solid ellipses), and in the case of IG dust with Rv = 9.5 (dashed ellipse). For clarity only the 68.3% level 
has been plotted for the dust distribution. 



y /dof xjjdof 



SN 2000fr 


no dust, (Q M , fi A )=(0.25,0.75) 
dustfl v = 9.5, (£2 M ,£2 A )=(1,0) 
dust/? v , = 4.5, (n M ,fi A )=(l,0) 


2.33/4 
4.24/4 
5.29/4 


2.08/3 
3.94/3 
4.50/3 


SN 1999ff 


no dust , (n M ,£l A )=(0.25,0.75) 
dustfl v = 9.5, (n M ,£2 A )=(l,0) 
dust/? v , = 4.5, (n M ,n A )=(i,0) 


6.05/5 
4.69/5 
4.31/5 


3.93/2 
2.06/2 
1.89/2 


SNe combined 


no dust , (n M ,fl A )=(0.25,0.75) 
dust2? v = 9.5, (£2 M ,£2 A )=(1,0) 
dustfl v = 4.5, (n M ,n A )=(l,0) 


8.38/9 
8.93/9 
9.59/9 


6.01/5 
6.00/5 
6.39/5 



Table 10 x 2 computed for the 3 different models and colours 
for each of the supernovae and for their combination. 

supernova-by-supernova or property -by-property basis, always 
finding counterparts of distant events in the local sample. 



In this work we compared the colours of nearby and dis- 
tant supernovae (primarily to test presence of "grey" dust). 
Although the size of the high redshift sample is very lim- 
ited, our results do not show any evidence for evolution in the 
colours of SNe la. Furthermore, the correlation found between 
the intensity of the secondary peak of /-band light curve and 
the supernova luminosity give an independent way of testing 
for evolution. The restframe /-band light curve of the high red- 
shift supernovae were all best fitted by templates showing a 
prominent second peak, i.e. inconsistent with the intrinsically 
underluminous supernovae. Note that the data presented here 
for SN 2000fr show for the first time a case where the sec- 
ondary peak is unambiguously evident in the data even prior 
to the light curve fit. Table lists the A^ 2 for the fit of the 
high redshift SNe to the templates of the two underluminous 
SN 1991bg and SN 1997cn, relative to the best fit. The x 2 val- 
ues are significantly larger than the best fit value. 
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n SN 1991bg SN 1997cn 
SN 2000fr 3 241)3 21/73 
SN 1999ff 2 3.28 2.72 



Table 1 1 Ajf 2 for the fit of the high redshift SNe to the templates 
of the two underluminous SNe relative to the best fits (which 
are "normal" SN templates), n is the number of data points used 
in the one-parameter fit (see discussion in section l4"3V 

8. Summary and conclusions 

In this work we have investigated the feasibility and utility 
of using restframe /-band observations for cosmological pur- 
poses. 

We have developed a five parameter light curve fitting pro- 
cedure which was applied successfully to 42 nearby Type la 
supernovae. The fitted light curves were used to build a set of 
templates which include a broad variety of shapes. We have 
found correlations between the fitted parameters, in particular 
between the time of the secondary peak and the B-band stretch, 
s B . Moreover, a width-luminosity relation was found between 
the peak /-band magnitude and the B- and /-band stretches (s B 
and sj). 

We built a restframe /-band Hubble diagram using 26 
nearby supernovae at redshifts 0.01 < z < 0.1, and measured an 
r.m.s. of 0.24 mag, smaller than the uncorrected dispersion cor- 
responding to restframe B-band. The width-luminosity relation 
was used to reduce the r.m.s. to 0.17 + 0.03 mag (including 
measurement errors), corresponding to an intrinsic dispersion 
of 0.13 mag. Differences between the three data samples are 
also discussed. 

/-band measurements of one new high redshift supernova 
plus published data of another were used to extend the Hubble 
diagram up to z ~ 0.5. The restframe /-band light curves of the 
z ~ 0.5 supernovae were fitted with templates that were built 
from the nearby SNe la, as the five parameter fit method could 
not be used for the poorly sampled high redshift light curves. 
The peak /-band magnitude of the high redshift SNe was com- 
pared to three different sets of cosmological parameters. The 
"concordance model" of the universe, Qa)=(0.25,0.75), 
is formally found in better agreement with the data than the 
other models at the >2cr level. However, the small sample size 
does not yet allow strong conclusions to be drawn. 

Alternative explanations for the observed dimming of su- 
pernova brightness, such as the presence of grey dust in the 
IG medium or evolutionary effects in the supernova properties 
have also been addressed. Both the /-band Hubble diagram and 
multi colour photometry have been used for testing grey dust. 
Although no firm limits on the presence of grey dust could be 
set, this study shows that with higher statistics, the restframe /- 
band measurements could provide useful information on cos- 
mological parameters, including tests for systematic effects. 
A Monte Carlo simulation indicates that a sample of at least 
20 well observed SNe la would be enough for setting limits 
through the multi-colour technique used in this paper. A simi- 
lar technique, using QSO instead of SNe la, was successfully 
used by Mortsell & GoobaJ d2003l) to rule out grey dust as be- 



ing the sole explanation for the apparent faintness of SNe la at 
z ~ 0.5. 

Possible systematic uncertainties affecting the restframe /- 
band Hubble diagram are discussed. Some sources are identi- 
fied, for instance the different methods applied for fitting the 
low and the high redshift samples, selection effects for bright 
objects due to the limiting magnitude of the search campaign, 
as well as uncertainties in the ^-correction calculations due to 
the presence of the Ca IR triplet feature in the near infrared 
region of the SN spectra. However, these systematic uncertain- 
ties differ from the ones that could affect the restframe B-band 
Hubble diagram. 

Restframe /-band observations of distant SNe la are feasi- 
ble, useful and complementary to the already well established 
observations in the B-band. 
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